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PARTI 


Ab  Initio  Quantum  Chemical  Design  of  Single  Supermolecule  Photoactive  Machines  and 

Molecular  Logical  Devices 

Arvydas  TAMULIS,  Jelena  TAMULIENE 

Institute  of  Theoretical  Physics  and  Astronomy,  A.  Gostauto  12,  2600  Vilnius,  Lithuania 

ABSTRACT 

It  were  designed  few  new  two  and  three  variable  molecular  logic  machines  which  possess 
excitation  energies  in  infra  red  (IR)  region,  as  well  as  selective  two  variable  molecular  logic 
devices  and  binary  two  and  three  variable  logic  machines. 

OR  logically  controlled  molecule  machine  designed  from  benzene-N=N-benzene-N02  with 
electron  acceptor  fragment  -(PI1-NO2)  (moving  part  after  excitation)  and  two  photoelectron  donor 
parts:  dithicno|3,2-b:2\3’-d]thiophenc  and  ferocene  molecules  show  some  not  intensive 
absorption  in  IR  region.  Applied  Density  Functional  Theory  (DFT)  -  Time  Dependent  (TD) 
method  and  our  modified  visualization  program  showed  small  charge  transfer  in  first  and  second  IR 
excited  stales. 

OR  logically  controlled  molecule  machine  designed  from  thiobenzene-N=N-benzene-N02  with 
electron  acceptor  fragment  -(Ph-N02)  (moving  part  after  excitation)  and  two  photoelectron  donor 
parts:  dithieno|3,2-b:2\3’-d]thiophene  and  ferocene  molecules  was  optimized  by  HF  method  and 
shows  characteristic  (similar  to  tetrahedral)  valence  angles  near  the  S  atom  inserted  in  the  six 
membered  ring. 

It  was  designed  three  variable  logically  controlled  molecule  machine  for  the  photoinduced 
moving  in  3Dimcnsional  zeolite  surfaces  and  construction  of  3D  addressing  memory. 

Design  and  optimization  of  geometry  by  HF/6-31G  was  performed  for  two  different  selective 
two  variable  molecular  logic  machines. 

It  was  designed  two  different  binary  two  and  three  variable  molecular  logic  machines  moving  in 
accordance  to  relaxation  energies  due  to  various  photo  excitations. 

Wc  performed  design  and  calculations  of  three  variable  molecular  logical  devices  based  on 
organic  electron  donor:  dilhicno[3,2-b:2\3’-d]thiophene,  tetrathiofulvalene  (TTF)  and  ferocene 
(C ] ()H  1  ()Fc)  and  electron  acceptor  molecules:  l,3-bis(dicyanomethylidenc)indanc  and  endohcdral 
fullcrenc  ErSc2N@C80  substituted  derivative  ErSc2N@(CH2CH2NH)C80>  electron  donor- 
bridgc-elcclron  acceptor  dyads  and  triads  including  electron  donor  and  acceptor  molecules  joined 
with  -CH=CH-  bridge.  Design  of  new  scries  molecular  implementations  of  two  variable  logic- 
functions:  AND  (NAND),  OR  (NOR)  is  based  on  quantum  chemical  ab  initio  Hartree-Fock 
( I  IF)/6— 3 1 1G  basis  geometry  optimization  procedure  which  shows  that  our  newly  designed  logical 
gates  based  on  ErSc2N@CgO  substituted  derivative  ErSc2N@(CH2CH2NH)C80  are  more  stable 
rather  than  based  on  ErSc2N@C80  substituted  derivative  ErSc2N@(CH2)C80- 

It  was  performed  design  and  calculation  of  biliverdine  derivatives  which  possess  from  one  to  six 
Quantum  Bits  generating  by  proton  NMR  for  closed  and  open  shell  species  (including  Fermi 
Contact  interaction  but  not  including  other  relativistic  effects).  Theory  of  our  obtained  different 
relativistic  terms  for  open  shell  paramagnetically  shifted  proton  NMR  calculations  and  formulae 
prepared  for  programming  and  implementation  to  the  DALTON  package  are  given  in  Appendix  1 
of  this  Report  04. 
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Set  consisting  from  six  various  ErSc2N@(CH2CH2NH)CH()  derivatives  where  H  is  replaced  by 
Cl:  ErSc2N@(CH2CH2NCl)C80,  ErSc2N@(CH2CClHNH)C80,  ErSc2N@(CClHCH2NH)C80, 
ErSe2N@(CH2CCl2NH)C80  etc.  were  designed  and  optimized  by  HF/Watanabe  in  order  to  get 
quantum  computing  elements  possessing  six  QuBits  based  on  six  different  EPR  signals. 

It  was  designed  molecular  switch  based  on  proton  transfer  in  molecular  complex  device 
composed  from  cuprate-mediate  oxidative  P-0  coupling  of  white  phosphorus  and  alcohols. 


Keywords:  quantum  mechanical  time-dependent  method,  charge  transfer  visualization,  infrared 
photoactive  molecular  logic  machines,  basic  elements  of  molecular  classical  and 
quantum  computers,  molecular  proton  switch 


Correspondence:  Email:  Tamulis@ITPA.lt;  WWW:  http://www.itpa.lt/~tamulis/;  Fax:  +370-2-225361 


1.  INTRODUCTION 

Analysis  of  presentations  of  Polish  International  Conference  "Towards  Molecular  Electronics" 
held  in  Srcm  (Poland),  25-30  June  2001  shows  that  in  previous  soviet  occupied  countries  rising 
new  generation  of  scientist  working  in  molecular  nano-technologies  based  mainly  on  organic 
nanoslruclurcd  compounds.  Our  student  Vykintas  Tamulis  presented  our  research  and  contacted 
with  most  interesting  speakers  of  thisTME’01  conference. 

Dr.  P.  Byszcwski  ct  al.  (Poland)  reported  concerning  molecular  memory  based  on  fullerenc  C60 

organomctallic  cycloaducls.  New  chemical  modifications  of  fullerene  C6()  adding  amines  presented 
Prof.  J.  J.  Danger  ct  al.  (Poland).  Dr.  A.  Lachinov  et  al.  (Russia)  working  with  Prof.  W.R.  Salaneck 
(Sweden)  presented  plenty  of  electronic  devices  including  LED  based  on  wide-gap  polymeric 
materials.  We  suggested  to  add  two  new  radicals  to  their  interesting  five  membered  rings 
monomers  (possessing  oxygen  atoms  that  allow  easy  remove  electrons)  two  new  radicals: 
photodonors  Dithienof3,2-b:2’,3’-d]thiophene  and  Ferocene  (ClOHlOFe).  We  planed  to  start 
calculate  these  monomers:  optimize  geometry  by  Hartree  Fock  6-31 1G**  basis  and  calculate 
spectra  and  charge  transfer  in  different  excited  states  by  Density  Functional  Theory  model 
B3PW9 1/6-31 1G**  using  Time  Dependent  method.  Vykintas  Tamulis  will  optimize  geometry  in 
all  excited  slates  by  Hartree  Fock  CIS  or  other  methods.  Later  he  will  add  to  these  five  membered 

rings  monomers  some  good  electron  acceptors  molecules:  TCNQ  ,  fullerenc  Ceo,  etc.  and  repeat  all 
mentioned  calculations  once  again  in  order  to  get  new  family  of  effective  molecular  devices. 

We  have  received  interesting  proposal  to  chair  SPIE  conference  "Phlochemical  Molecular  & 
Supramolccular  Devices"  to  be  hold  in  2002  17-21  June  Messe  Frankfurt.  It  is  foreseen  lecture  of 
J.-M.  Lchn  and  several  sessions:  Photonics  signals,  photocleavable  cryptates,  molecular 
switching/EFISH  (electric  field  induced),  molecular  wires,  optical  memory,  elementary  processes 


in  molecular  computing,  etc. 

Our  research  presented  in  this  Report  04  coincide  with  themes  of  above  mentioned  conferences. 


Mainly  we  arc  performing  design  of  photoinduced  electron  transfer  molecular  devices. 


Electron  transfer  in  molecular  devices  constructed  from  electron  donors,  electron  acceptors  and 
insulating  bridges  molecular  fragments  can  be  divided  in  two  different  processes: 

/ )  direct  electron  transfer ,  then  electron  excited  by  radiation  from  molecular  orbital  localized 
on  donor  fragment  of  molecule  to  molecular  orbital  localized  on  electron  acceptor  fragment 
2)  indirect  electron  transfer ,  then  electron  excited  from  electron  donor  molecular  orbital  to 
empty  same  electron  donor  molecular  orbital,  and  electron  transfer  occurs  as  non-radiate 
relaxation  of  electron  donor  excited  state 

Both  processes  can  by  described  using  Marcus  electron  transfer  theory.  This  theory  shows  that 
electron  transfer  rate: 


v  =  —\V?(FC) 
h  M  (1) 

where  V  -  electron-electron  coupling  matrix  element,  (FC)  -  Frank-Condon  facto.  Franc-Condon 
factor  in  first  approximation: 


_Ea 

(FC)  =  c~kT 


E,  = 


(AG  +  Xf 
41 


where  X  -  reorganization  energy  of  molecular  structure.  In  case  of  adiabatic  direct  electron  transfer 
we  have: 


Excited  state 


From  this  picture  sec  how  to  obtain  parameters  for  (FC)  approximation.  Electron-electron 
coupling  matrix  elements  can  be  calculated  using  potential  energy  surface  search  or  applying 
external  electric  field  to  bring  donor  and  acceptor  levels  into  resonance.  Currently  we  are 
performing  calculations  of  direct  electron  transfer  rates  using  external  electric  field  method, 
because  it  requires  less  computational  efforts. 
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2,  RESULTS  AND  DISCUSSIONS 


2.1.  Design  of  light  induced  two  and  three  variable  logically  controlled  molecular  machines, 
3D  addressing  memory,  selective  logic  molecular  machines  and  binary  logically  controlled 

molecular  machines 


Two  and  three  variable  OR  logically  controlled  molecular  machines  are  composed  from 
optimized  by  B3PW9 1/6-3 1 1G**  separate  azo-dye  molecules:  benzene-N=N-benzene-N02, 
lhiobcnzcnc-N=N-bcn/.cnc-N02  and  photodonor  molecules:  dithieno[3,2-b:2\3’-d]thiophcnc, 
fcroccnc  (CjqHioFc)  (formula  is  C12H9N),  tetrathiofulvalene  (TTF),  carbazole  (Cz)  and 
phcnylcncdiaminc  (PhD A).  These  molecular  machines  were  designed  based  on  found  internal 
molecular  motions  of  azo-dye  D03  and  azo-benzene  molecules  [1-4]  (see  also  our  EOARD 
Reports  01-03).  Photodonor  molecules  (thiophene,  ferocene,  TTF,  Cz,  PhD  A)  were  attached  to 
azo-dye  molecules:  benzenc-N=N-benzene-NC>2,  thiobenzene-N=N-benzene-NC>2  by  using 
-HC=CH-  bridges  (sec  Figure  s  1,  7,  8,  12,  13). 

2.1.1.  Design  of  light  induced  two  variable  logically  controlled  molecular  machines  based  on 

ferocene,  thiophene  and  D03 

OR  molecular  device  composed  from  good  photo  electron  donors  ferocene  and  thiophene  (see 
Figure  1)  should  shift  excitation  energies  to  red  or  even  IR  region.  In  the  case  if  the  thiophene  or 
fcroccnc  molecules  arc  excited  by  light,  electron  should  jump  passing  the  -N=N-  bridge  will 
change  the  geometry  of  the  supermolecule.  After  the  dissipation  of  excited  energy  on  the  surface 
the  supermolecule  will  return  to  its  ground  state  geometry  moving  on  the  surface.  Movement  of  a 
random-walker  it  is  possible  if  thiophene  or  ferocene  fragments  are  excited  therefore  this  device 
should  be  two  variable  logically  OR  controlled  molecular  random-walker. 


Figure  1 .  Two  variable  OR  logically  controlled  molecule  machine  designed  from  benzene-N=N- 
bcn/,cnc-N02  with  electron  acceptor  fragment  -(PI1-NO2)  (moving  part  after  excitation)  and  two 
pholoclcclron  donor  parts:  dithieno[3,2-b:2’,3’-d]thiophene  and  ferocene  molecules. 

The  spectrum  of  this  two  variable  OR  logically  controlled  single  supermolecule  device  was 
calculated  by  DFT-B3PW91  TD  using  6-31 1G  basis  set  (see  Table  1). 

Table  1 .  Excitation  energies  and  oscillator  strengths  of  two  variable  OR  device  calculated  by  DFT 
B3PW91-TD  6-31 1G  basis  set. 

Excitation  energies  and  oscillator  strengths: 

Excited  State  1:  Singlet  0.8760  eV  1415.30  nm  f=0.0002 
168  ->172  .18760 

168  ->173  .25566 

168 ->174  -.27833 

168 ->175  .21666 

168 - >177  .18935 

169  ->172  .11177 

169  ->173  .14934 

169 - >175  -.42337 

169  ->177  .11911 

Excited  State  2:  Singlet  0.8774  cV  1413.14  nm  f=0.0005 

168  ->172  -.11134 

168 ->173  -.14833 

168 -> 1 75  .42167 

168 - >177  -.11709 

169 - >172  .18419 

169  ->173  .25290 

169  ->174  -.27837 
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169  ->175 

.22201 

169  ->177 

.19218 

Excited  State 

3:  Singlet 

1.3294  eV 

932.61  nm 

f=0.0006 

163  ->172 

.13754 

163  ->173 

.19337 

163  ->174 

-.18590 

163  ->175 

-.10026 

163  ->177 

.15819 

168  ->175 

-.35464 

169  ->172 

.12293 

169  ->173 

.15703 

169 ->174 

-.15519 

169 ->177 

.10524 

Excited  State 

4:  Singlet 

1.3380  eV 

926.66  nm 

f=0.0001 

163  ->175 

-.32656 

168  ->172 

.13235 

168  — > 1 7  3 

.16668 

168  ->174 

-.16094 

168  ->177 

.10569 

169 ->175 

.33869 

Excited  State 

5:  Singlet 

1.7353  eV 

714.49  nm 

f=0.0014 

169  ->170 

.70408 

Excited  State 

6:  Singlet 

1.7874  eV 

693.66  nm 

f=0.0000 

168  ->170 

.70576 

Excited  State 

7:  Singlet 

2.0972  eV 

591.18  nm 

f=0.0002 

164  ->170 

.60450 

164  ->171 

.16532 

165  ->170 

.17304 

Excited  State 

8:  Singlet 

2.2812  eV 

543.51  nm 

f=0,0000 

163  ->172 

.22870 

163  ->173 

.30490 

163  ->174 

-.30564 

163  ->177 

.22008 

168  ->175 

.16429 

Excited  State 

9:  Singlet 

2.3343  eV 

531.15  nm 

f=0.0000 

163  ->175 

-.54793 

168  ->172 

-.10295 

169  ->175 

-.16357 

Excited  State 

10:  Singlet 

2.3925  eV 

'  518.22  nm 

f=0.0083 

167  ->170 

.68181 

Excited  Stale 

11:  Singlet 

2.5855  eV  479.54  nm 

f=0.0391 

166 ->170 

.68680 

167  ->170 

.10215 

Excited  Slate 

12:  Singlet 

2.8212  eV 

439.47  nm 

f=0.0006 

169 ->171 

.70219 

Excited  State 

13:  Singlet 

2.8722  e\ 

r  431.67  nm 

i  f=0.0000 

168  ->171 

.70216 

Excited  State 

14:  Singlet 

2.9414  e\ 

'  421.50  nm 

i  f=0.0127 

164  ->170 

-.19180 

165  ->170 

.67276 

Excited  State 

15:  Singlet 

3.0164  eV 

411.02  nm 

f=0.0011 

163  ->170 

.70524 

From  Tabic  1  il  possible  to  see  that  this  OR  device  possesses  some  not  intensive  absorption  lines 
in  infrared  (IR)  region. 

Our  calculated  B3PW91-TD/6-31 1G  charge  in  ground  (see  Figure  2)  and  charge  transfer  in 
first  and  second  excited  states  (sec  Figures  4,  6)  using  calculation  files  are  visualized  by  our 
modified  visualization  and  accounting  software  given  in  Appendix  1  of  our  Report  3.  Black  circles 
correspond  to  electron  removing  values  and  white  circles  correspond  to  electron  income  values. 
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Figure  2.  Charge  distribution  in  two  variable  OR  logical  device  in  ground  state. 


Figure  4.  Charge  transfer  i 
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two  variable  OR  logical  device  in  first  excited  state. 


Figure  5.  Charge  distribution  in  two  variable  OR  logical  device  in  second  excited  state. 


Figure  6.  Charge  transfer  in  two  variable  OR  logical  device  in  second  excited  state. 

Total  charge  transfer  in  first  and  second  excited  states  are  small:  equal  to  0.174  and  0.173 
electron  charge  because  the  oscillator  strengths  of  these  two  IR  transitions  are  weak:  f=0.0002  and 
1=0.0005  (sec  Table  1). 

The  further  research  in  this  field  is  continuingB3PW91-TD  calculations  in  other  excited  states 
up  it)  near  UV  region  searching  the  charge  transfer  from  two  different  inputs  of  this  single 
supcrmolcculc  OR  molecular  logic  gate-molecular  machine. 

Another  two  variable  OR  molecular  device  composed  from  good  photo  electron  donors  feroeene 
and  thiophene  and  inserted  S  atom  in  the  benzene  ring  (see  Figure  7)  should  significantly  shift 
excitation  energies  to  red  or  IR  region..  In  the  case  if  the  thiophene  or  feroeene  molecules  are 
excited  by  light,  electron  should  jump  passing  the  -N=N-  bridge  will  change  the  geometry  of  the 
supcrmolcculc.  After  the  dissipation  of  excited  energy  on  the  surface  the  supermolecule  will  return 
to  its  ground  state  geometry  moving  on  the  surface.  Movement  of  a  random-walker  it  is  possible  if 
thiophene  or  feroeene  fragments  are  excited  therefore  this  device  should  be  two  variable  logically 
OR  controlled  molecular  random-walker. 


Figure  7.  Two  variable  OR  logically  controlled  molecule  machine  designed  from  thiobenzenc- 
N=N-bcnzcnc-N02  with  electron  acceptor  fragment  -(PI1-NO2)  (moving  part  after 
excitation)  and  two  photoclcctron  donor  parts:  dithieno[3,2-b:2’,3’-d]thiophene  and 
feroccnc  molecules 

The  geometry  of  two  variable  OR  device  presented  in  Figure  7  was  optimized  by  HF  using  SVP 
basis  and  gave  characteristic  close  to  tetrahedral  dihedral  angles  (104,  100,  108)  around  the  inserted 
to  benzene  ring  S  atom  (the  similar  angles  we  obtained  in  CdS  clusters  (sec  our  EOARD  Reports 
02  and  03).  This  is  because  of  strong  interaction  of  lone  pairs  of  electrons  of  S  atom. 

DFT-TD  B3PW91  with  SVP  basis  calculations  of  electronic  structure  of  this  OR  device  shows 
stability  even  near  S-bcnzcnc  ring  atoms.  Overlapping  population  (OP)  positive  values  for  bonds 
-C-S-C-  arc  equal  to  0.405  and  0.445  and  OP  for  bond  =S-N=  is  equal  to  0.765. 

Further  development  of  this  research  is  calculations  of  spectra  and  charge  transfer  using 
B3PW9 1/6-31 1G-TD  method. 

2.1.2.  Design  of  light  induced  three  variable  logically  controlled  molecular  machines  and  3D 

addressing  memory 

Three  variable  OR  molecular  device  composed  from  good  photo  electron  donors  TTF,  ferocenc 
and  thiophene  (sec  Figure  8)  should  allow  to  move  this  photo  driven  charge  transfer  molecular 
machine  in  three  dimensional  surfaces  (for  example,  in  zeolites  (see  Figure  9)).  In  the  case  if  the 
TTF.  thiophene  or  feroccnc  molecules  arc  excited  by  light,  electron  should  jump  passing  the 
-N=N-  bridge  will  change  the  geometry  of  the  supermolecule.  After  the  dissipation  of  excited 
energy  on  the  zeolite  surface  the  supermolecule  will  return  to  its  ground  state  geometry  moving  on 


the  zeolite  surface.  Movement  of  a  random-walker  it  is  possible  if  TTF,  thiophene  or  ferocenc 
fragments  arc  excited  therefore  this  device  should  be  two  variable  logically  OR  controlled 
molecular  random-walker. 


Figure  8.  Three  variable  OR  logically  controlled  molecule  machine  designed  from  bcnzenc-N=N- 
ben/.cnc-N02  with  electron  acceptor  fragment  -(PI1-NO2)  (moving  part  after  excitation)  and  three 
pholoclcclron  donor  parts:  dithieno[3,2-b:2\3’-d]thiophene,  ferocene  and  TTF  molecules. 

Movement  of  this  three  variable  OR  molecular  logic  device  it  is  possible  to  control  exciting 
TTF,  thiophene  or  ferocenc  photoelectron  donor  molecules  by  certain  order  which  allow  to  move 
OR  device  to  certain  space  in  zeolite  (see  Figure  9)  that  should  result  the  fixation  of  signal  in  three 
dimensional  addressing  memory. 


Figure  9.  Zeolite-like  three  dimensional(3D)  addressing  memory  where  moves  three  variable  OR 
logieally  controlled  molecule  machine  designed  from  benzene-N=N-benzene-N02  with  electron 
acceptor  fragment  -(Ph-N02)  (moving  part  after  excitation)  and  three  photoelectron  donor  parts: 
dilhicno|3,2~b:2\3’-d|thiophcnc,  fcrocenc  and  TTF  molecules. 

2.1.3.  Design  of  light  induced  selective  logic  molecular  machines 

Depending  on  excitation  energy  charge  transfer  might  be  to  the  different  branches  of  selective 
two  variable  OR  logically  controlled  molecule  machine  (see  Figure  10)  designed  from  two 
photoclcctron  donors:  bcn/cnc-CH=CH-bcnzcne-NH2,  TTF  and  two  electron  acceptor  molecules: 
-bcn/,cnc-CH=CH-TCNB,  -benzene-N=N-benzene-N02  •  Electron  acceptor  fragment  -(Ph- 
NO2)  is  moving  part  after  excitation.  Therefore  we  should  have  two  combinations  of  two  variable 
OR  molecular  logic  machines  on  the  same  supermolecule  moving  by  two  different  ways  depending 
on  certain  excitation  energies. 

Geometry  of  this  OR  molecular  was  optimized  by  HF  using  6-3 1G  basis. 
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Figure  10.  Selective  two  variable  OR  logically  controlled  molecule  machine  designed  from  two 
photoclcctron  donors:  bcnzcnc-CH=CH-benzene-NH2,  TTF  and  two  electron  acceptor  molecules: 
-bcn/.cnc-CH=CH-TCNB,  -benzene-N=N-benzene-N02  •  Electron  acceptor  fragment  -(Ph- 
NO2)  is  moving  part  after  excitation. 


Another  kind  of  selective  two  variable  OR  logically  controlled  molecule  machine  designed  from 
two  better  photoclcctron  donors:  Fcroccnc,  TTF  and  two  electron  acceptor  molecules:  1,3- 
bis(dicyanoinethylidcnc)indanc,  -bcnzcnc-N=N-bcnzene-N02  .  Electron  acceptor  fragment 
-(PI1-NO2)  is  moving  part  after  excitation.  We  expect  to  excite  this  selective  two  variable 
molecular  logic  device  by  red  or  even  IR  energies. 


Figure  1 1 .  Selective  two  variable  OR  logically  controlled  molecule  machine  designed  from  two 
photoclcetron  donors:  Ferocene,  TTF  and  two  electron  acceptor  molecules: 

1 ,3-bis(dicyanomcthylidcnc)indanc,  -bcnzene-N=N-benzene-N02  .  Electron  acceptor  fragment 
-(PI1-NO2)  is  moving  part  after  excitation. 


The  further  development  of  these  both  selective  OR  molecular  devices  is  performing  calculating 
spectra  and  charge  transfer  traces  using  B3PW9 1/6-3 1 1G  -TD  method. 


2.1,4.  Design  of  light  induced  two  and  three  variable  binary  logically  controlled  molecular 

machines 

Supramolcculcs  composed  from  two  and  three  OR  logically  controlled  molecular  machines 
should  move  in  accordance  of  their  exiting  spectra  corrected  by  Van  dcr  Waals  interactions.  For  the 
evaluation  of  charge  transfer  reorganization  energies  firstly  was  calculated  spectrum  and  charge 
redistribution  of  OR  logical  device  composed  from  Carbazole  attached  to  D03  molecule  (see 
Figure  12). 

Geometry  optimization  and  spectrum  was  calculated  by  AMI  method  (see  Figure  12)  and  Table 


Figure  12.  Two  variable  OR  logically  controlled  molecule  machine  designed  from  two 
pholoclcctron  donors:  Carbazole  and  electrondonor  fragment  -benzene-NHo  and  moving 
clcctronacccptor  part — N=N-bcnzcnc-NC>2. 


l’uble  2.  Several  details  on  in  Figure  12  placed  OR  logical  device  calculated  by  AMI 

STATE  E(EV)  E(CM-l)  L(NM')  OSC  TX  TY  TZ  TM 
I  2.585  20847.6  479.7  0.00002  -0.04  0.00-0.01  0.04 

Change  of  ihc  total  charge  on  FRAGMENT  1  (atoms  1-  14  )  0.199 
Change  of  the  total  charge  on  FRAGMENT  2  (atoms  15-28  )  -0.179 

Change  of  the  total  charge  on  FRAGMENT  3  (atoms  30-  50 )  -0.020 

Dipole  moment  X  Y  Z  Total 
Ground  state  -8.495  0.422  1.222  8.593 

Change  by  exc.  -9.526  1.409  -0.057  9.630 

Excited  state  -18.021  1.831  1.165  18.151 


STATE  E(EV)  E(CM-l)  L(NM)  OSC  TX  TY  TZ  TM 
2  3.340  26941.0  371.2  0.85740  -8.22  0.21  0.28  8.23 

Change  of  the  total  charge  on  FRAGMENT  1  (atoms  1-  14  )  0.156 
Change  of  the  total  charge  on  FRAGMENT  2  (atoms  15-28  )  -0.583 
Change  of  the  total  charge  on  FRAGMENT  3  (atoms  30-  50  )  0.384 

Dipole  moment  X  Y  Z  Total 
Ground  state  -8.495  0.422  1.222  8.593 

Change  hv  exc.  -28.253  -10.436  1.603  30.161 

Excited  state  -36.748  -10.014  2.825  38.192 


STATE  E(EV)  E(CM-l)  L(NM)  OSC  TX  TY  TZ  TM 
3  3.669  29590.3  337.9  0.35797  2.38  4.42-0.74  5.07 

Change  of  the  total  charge  on  FRAGMENT  1  (atoms  1-  14  )  0.129 
Change  of  the  total  charge  on  FRAGMENT  2  (atoms  15-  28  )  -0.085 
Change  of  the  total  charge  on  FRAGMENT  3  (atoms  30-  50  )  -0.036 
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Dipole  moment  X  Y  Z  Total 

Ground  state  -8.495  0.422  1.222  8.593 

Change  by  cxc.  -3.630  1.941  -0.452  4.141 

Excited  state  -12.125  2.363  0.770  12.377 


STATE  E(EV)  E(CM-l)  L(NM)  OSC  TX  TY  TZ  TM 
4  3.741  30171.0  331.4  0.40258  -5.16  1.33  -0.06  5.33 

Change  of  the  total  charge  on  FRAGMENT  1  (atoms  1-  14  )  0.233 
Change  of  the  total  charge  on  FRAGMENT  2  (atoms  15-  28  )  -0.430 
Change  of  the  total  charge  on  FRAGMENT  3  (atoms  30-  50 )  0.171 

Dipole  moment  X  Y  Z  Total 
Ground  state  -8.495'  0.422  1.222  8.593 

Change  by  cxc.  -20.169  -4.098  0.850  20.598 

Excited  state  -28.664  -3.675  2.072  28.973 


STATE  E(EV)  E(CM-l)  L(NM')  OSC  TX  TY  TZ  TM 

5  4.033  32527.0  307.4  0.39450  -1.80-4.70  0.66  5.08 

Change  of  the  total  charge  on  FRAGMENT  1  (atoms  1-  14  )  0.051 
Change  of  the  total  charge  on  FRAGMENT  2  (atoms  15-  28  )  -0.089 
Change  of  the  total  charge  on  FRAGMENT  3  (atoms  30-  50  )  0.053 

Dipole  moment  X  Y  Z  Total 
Ground  state  -8.495  0.422  1.222  8.593 

Change  by  exc.  -5.029  -0.731  -0.158  5.085 

Excited  state  -13.524  -0.309  1.064  13.569 


STATE  E(EV)  E(CM-l')  L(NM)  OSC  TX  TY  TZ  TM 
6  4.154  33504.3  298.5  0.03702  0.13  1.53  -0.08  1.53 

Change  of  the  total  charge  on  FRAGMENT  1  (atoms  1-  14  )  -0.066 
Change  of  the  total  charge  on  FRAGMENT  2  (atoms  15-  28  )  0.024 
Change  of  the  total  charge  on  FRAGMENT  3  (atoms  30-  50  )  0.037 

Dipole  moment  X  Y  Z  Total 
Ground  state  -8.495  0.422  1.222  8.593 

Chanue  by  cxc.  -2.758  -1.213  0.184  3.018 

Excited  state  -11.253  -0.791  1.406  11.368 


STATE  E(EV)  E(CM-l)  L(NM)  OSC  TX  TY  TZ  TM 
7  4.272  34458.7  290.2  0.07989  0.76-0.13  2.08  2.22 

Change  of  the  total  charge  on  FRAGMENT  1  (atoms  1-  14  )  -0.083 
Change  of  the  total  charge  on  FRAGMENT  2  (atoms  1 5-  28  )  -0.01 1 
Change  of  the  total  charge  on  FRAGMENT  3  (atoms  30-  50  )  0.135 

Dipole  moment  X  Y  Z  Total 
Ground  state  -8.495  0.422  1.222  8.593 

Change  by  exc.  -2.062  -2.939  0.230  3.598 

Excited  state  -10.557  -2.517  1.452  10.949 


STATE  E(EV)  E(CM-l)  L(NM)  OSC  TX  TY  TZ  TM 

8  4.343  35029.3  285.5  0.15366  0.86-2.93  0.18  3.05 

Change  of  the  total  charge  on  FRAGMENT  1  (atoms  1-  14  )  0.137 
Change  of  the  total  charge  on  FRAGMENT  2  (atoms  15-  28  )  -0.325 
Change  of  the  total  charge  on  FRAGMENT  3  (atoms  30-50)  0.168 

Dipole  moment  X  Y  Z  Total 
Ground  state  -8.495  0.422  1.222  8.593 

Chantre  by  cxc.  -15.784  -4.301  0.688  16.373 

Excited  state  -24.278  -3.879  1.910  24.660 

STATE  E(EV)  E(CM-l)  L(NM)  OSC  TX  TY  TZ  TM 

9  4.413  35589.4  281.0  0.04703  0.54-0.46  1.52  1.68 


Change  of  the  lolal  charge  on  FRAGMENT  1  (atoms  1-  14  )  -0.017 
Change  of  the  total  charge  on  FRAGMENT  2  (atoms  15-  28  )  -0.025 
Change  of  the  total  charge  on  FRAGMENT  3  (atoms  30-  50  )  0.054 

Dipole  moment  X  Y  Z  Total 
Ground  state  -8.495  0.422  1.222  8.593 

Change  by  exc.  -1.911  -1.267  -0.097  2.294 

Excited  state  -10.405  -0.845  1.125  10.500 


STATE  E('EV)  E(CM-l)  L(NM)  OSC  TX  TY  TZ  TM 
10  4.624  37297.7  268.1  0.02131  0.41  1.00-0.24  1.10 

Change  of  the  total  charge  on  FRAGMENT  1  (atoms  1-  14 )  0.145 

Change  of  the  total  charge  on  FRAGMENT  2  (atoms  15-  28  )  -0.413 

Change  of  the  total  charge  on  FRAGMENT  3  (atoms  30-  50  )  0.245 

Dipole  moment  X  Y  Z  Total 
Ground  state  -8.495  0.422  1.222  8.593 

Chanuc  by  exc.  -21.779  -6.344  1.016  22.707 

Excited  state  -30.274  -5.922  2.238  30.929 


Figure  13.  Binary  supramolecule  composed  from  two  variable  OR  logically  controlled  molecule 
machine  designed  from  two  photoelectron  donors:  Carbazole  and  electrondonor  fragment 


-bctr/.cnc-NH2  and  moving  clectronacceplor  part — N=N-benzene-N02  and  parallely  displaced 
from  three  variable  OR  logically  controlled  molecule  machine  designed  from  three  photoelectron 
donors:  Carbazole,  PhDA  and  electrondonor  fragment  -benzene-NH2  and  moving  electronacceptor 
part — N=N-bcnzcnc-N02. 


For  the  evaluation  of  charge  transfer  reorganization  energies  it  were  designed  two  different 
supramolcculcs  composed  from  two  and  three  variable  OR  molecular  devices  (sec  Figures  13  and 
14). 
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Figure  14.  Binary  supramolecule  composed  from  two  variable  OR  logically  controlled  molecule 
machine  designed  from  two  pholoclectron  donors:  Carbazole  and  electrondonor  fragment 

-bcn/.cnc-NH2  and  moving  elcctronacceptor  part — N=N-benzene-N02  and  from  linearly 
displaced  three  variable  OR  logically  controlled  molecule  machine  designed  from  three 
photoclcctron  donors:  Carbazole,  PhDA  and  electrondonor  fragment  -benzene-NH2  and  moving 
elcctronacceptor  part — N=N-bcnzenc-N02. 


2.2.  Design  of  fullerene  derivatives  based  three  variable  molecular  logic  devices 


Wc  performed  design  and  calculations  of  molecular  logical  devices  based  on  organic  electron 
donors:  dithicno[3,2-b:2’,3’-d]thiophene,  ferocene  (CioHioFe),  and  tetrathiofulvalene  (TTF)  and 
electron  acceptor  molecules:  1 ,3-bis(dicyanomethylidene)indane  and  endohedral  fullerene 
ErSc2N@CB()  substituted  derivative  ErSc2N@(CH2CH2NH)C80,  electron  donor-bridge-electron 
acceptor  dyads  and  triads  including  electron  donor  and  acceptor  molecules  joined  with  -CH=CH- 
bridge.  Design  of  new  series  molecular  implementations  of  three  variable  logic  functions:  AND 
(NAND),  OR  (NOR)  is  based  on  quantum  chemical  ab  initio  HF/6— 3 1 1 G  geometry  optimization 
procedure  which  shows  that  our  newly  designed  logical  gates  based  on  ErSc2N@C80  substituted 
dcrivaliveErSc2N@(CH2CH2NH)C8()  are  more  stable  rather  than  based  on  ErSc2N@C80 
substituted  derivative  ErSc2N@(CH2)C80- 

The  cluster  ErSc2N  that  exists  in  the  endohedral  fullerene  ErSc2N@C8  0  shows  magnetic 
activity  calculated  by  Restricted  Open  Shell  HF  method  using  Watanabe  (WTBS)  [5,  6]  basis  set  in 
doublet  electronic  state. 

In  Figure  15  it  is  presented  three  variable  Conjunction  (AND)  or  NAND  gates  designed  from 
three  electron  donor  fragments:  dithieno[3,2-b:2’,3’-d]thiophene,  ferocene  (CioHioFe),  TTF  and 
two  electron  acceptor  molecules:  ],3-bis(dicyanomethylidene)indane  and  fullerene 
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-CH2CH2NHC8O  derivative)  connected  with  -HOCH-  bridges.  dithieno[3,2-b:2\3’- 
d] thiophene,  ferocenc  and  TTF  molecules  possess  different  excitation  wavelengths  therefore  might 
be  three  different  inputs  to  the  AND  gate.  Let  us  imagine  that  we  are  measuring  the  additional 
electron  charge  on  fullerene  -CHCH2NHC8O  fragment.  The  present  state  of  experimental  physics 
allows  measuring  the  single  electron  charge  on  molecule  [7-10].  In  the  case  if  we  do  not  excite  no 
one  of  two  pholoclcelron  donor  molecules  on  the  input  of  this  gate  we  should  not  measure  the 
additional  electron  charge  on  the  output  -  fullerene  acceptor  fragment  -CHCH2NHC80-  Therefore 
in  this  case  we  should  have  response  on  the  output:  0. 


Figure  15.  Three  variable  AND  or  NAND  gate  designed  from  three  electron  donor  fragments 
and  two  electron  acceptor  molecules. 

In  the  case  if  we  should  excite  by  light  one  of  three  photo  electron  donor  molecules  we  also 
should  not  measure  additional  electron  charge  on  the  output  of  this  gate  -  fullerene 
-CH2CH2NHC8O  fragment  (because  to  our  quantum  chemical  calculations  cyano-indane  molecule 
is  good  acceptor  molecule  and  electron  hopping  from  the  cyano-indane  to  fullerene  should  not 
exist  or  will  be  with  small  probability).  Therefore  in  all  two  these  cases  we  should  have  responses 
on  the  output:  0.  In  the  ease  if  we  should  excite  two  photo  electron  donor  molecules  (two  inputs)  by 
all  possible  combinations:  1)  Thiophene  and  Ferocene;  2)  Thiophene  and  TTF;  3)  Ferocene  and 
TTF  we  should  measure  electron  on  the  output  -  fullerene  molecule  with  large  probability. 
Therefore  we  should  have  the  response  on  the  output:  1.  Summarizing,  in  all  four  different  cases  of 
this  two  variable  logic  gate  we  should  have  sequence  of  responses  on  output:  0,0, 0,1.  It  means  that 
this  our  Ml  logic  function  is  three  variable  conjunction  (AND). 
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Wc  should  have  MI  logic  function  NAND  in  the  same  device  in  the  case  if  we  excite 
dilhicnol  3,2— b:2’ ,3’— d Jthiophcne  and  fcrocene  molecules  on  the  input  and  measure  on  the  output 
the  spectrum  of  the  neutral  (without  additional  electron  charge)  fullerene  -CH2CH2NHC6O 
molecule  because  wc  should  have  the  sequence  on  output:  1,1, 1,0. 

In  Figure  16  it  is  presented  Disjunction  (OR)  or  NOR  gate  designed  from  three  electron  donor 
fragments:  dithieno|3,2-b:2\3’-d]thiophene,  ferocene,  TTF  and  fullerene  derivative 
-CH2CH2NHC8O  acceptor  fragment  connected  by  -HOCH-  bridges. 

Let  us  imagine  that  wc  do  not  excite  by  light  no  one  photoelectron  donor  molecule  (not  send 
signals  to  the  input)  of  this  MI  gate.  Then  we  do  not  measure  the  additional  electron  charge  on  the 
fullerene  fragment  -CHCH2NHC8O  (on  the  output  of  this  gate).  In  this  case  we  should  have 
response:  0.  In  the  case  if  wc  by  light  one  of  photoelectron  donor  molecules  (one  of  three  inputs: 
fcrocene.  TTF  or  thiophene)  electron  should  jump  from  photoelectron  donor  to  fullerene.  In  this 
case  wc  should  have  the  response:  1.  In  the  case  if  we  excite  both  photoelectron  donor  molecules  in 
all  combinations:  ferocene  or/and  thiophene;  ferocene  or/and  TTF,  thiophene  or/and  TTF  we 
should  have  response  on  fullerene  molecule:  1.  Summarizing  in  all  four  possible  cases  wc  should 
have  sequence  of  responses  on  the  output  fullerene  molecule:  0,1, 1,1.  It  means  that  this  our  MI 
logic  gate  is  disjunction  (OR). 

Wc  should  have  MI  logic  function  NOR  in  the  same  device  in  the  case  if  we  excite  ferocene 
or/and  thiophene;  fcrocene  or/and  TTF,  thiophene  or/and  TTF  molecules  on  the  inputs  of  the  same 
gale  but  measure  the  output  spectrum  of  the  neutral  (without  the  additional  electron  charge) 
fullerene  fragment  -CHCH2NHC8O  because  we  should  have  the  sequence  of  responses  on  the 
output  molecule:  1,0, 0,0. 
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Figure  16.  Three  variable  OR  or  NOR  gate  designed  from  three  electron  donor  fragments:  ferocene 
and  thiophene  and  fullerene  acceptor  fragment  -CHCH2NHC80- 


The  further  development  of  this  research  is  calculating  charge  transfer  in  these  three  variable 
molecular  logic  gates  due  to  different  excitation  energies.  We  expect  to  get  interesting  results  in 
reorganisation  of  geometry  of  inside  fullerene  clusters  due  to  different  charge  transfer 
reorganisation  energies. 
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2.3.  Design  of  Elements  of  Quantum  Computers  Based  on  Proton  NMR  and  on  EPR 
Spectra 

2.3. 1.  Electron  structure  and  proton  NMR  calculations  of  biliverdin  derivatives. 


Bi  liver  din  based  'H  NMR  6  QuBits  device  presented  in  Report  03  was  additionally  optimized 

tightly  using  HF  method  6-31 1G  basis  (sec  Figure  17)  in  order  to  get  as  much  possible  exact  C2 
symmetry  and  was  additionally  applied  Localized  Orbitals  analysis. 


2  3 


Figure  17.  Geometry  of  additionally  optimized  geometry  of  6  QuBits  quantum  computing  device 
very  close  to  C2  symmetry. 


Table  3.  Charges  on  atoms  calculated  by  HF  method  using  6-31 1G  basis: 
Total  atomic  charges: 


1  Cu 

0.931 197 

2  N 

-0.670304 

3  C 

0.606696 

4  N 

-0.721585 

5  0 

-0.401348 

6  C 

-0.277643 

7  C 

-0.105742 

8  C 

0.143673 

9  C 

0.032919 

10  C 

0.180590 

11  c 

-0.265247 

12  C 

-0.194847 

13  C 

0.128443 

14  C 

0.079271 

15  H 

0.208873 

16  H 

0.200848 

17  H 

0.176999 

18  11 

0.180520 

19  H 

0.181961 

20  H 

0.179866 

21  N 

-0.670304 

22  C 

0.606721 

23  N 

-0.721591 

24  0 

-0.401373 

25  C 

-0.277647 

26  C 

-0.105736 

27  C 

0.143704 

28  C 

0.032882 

29  C 

0.180628 

30  C 

-0.265268 

31  C 

-0.194818 

32  C 

0.128437 

33  H 

0.180526 

34  H 

0.208880. 

35  H 

0.200853 

36  H 

0.181961 

37  H 

0.177003 

It  is  possible  to  recognize  the  C2  symmetry  from  the  values  of  Table  3  having  atom  numciation 
in  Figure  17. 

In  order  to  ensure  which  electronic  state  of  6  Bits  quantum  computing  device  is  more  preferable 
it  Was  done  systematic  calculation  s  of  neutral  and  ionic  molecules  possessing  different  multiplicity. 
The  results  of  these  investigations  are  summarized  in  Tables  4  and  5. 


Table  4.  Total  energy  of  neutral  6  QuBits  device  possessing  C2  symmetry  in  various 


ilium  piitiiiu') 

Method 

Multiplicity 

Before  Opimization 

After  Optimization 

ROHF 

Singlet 

-2731.86 

-2731.86 

ROHF 

Triplet 

-2731.81 

-2731.82 

ROHF 

Pentet 

-2731.86 

-2731.89 

ROHF 

Septet 

-2731.73 

Method 

Multiplicity 

Energy 

B3PW9 1/6-3 1 1G 

1 

-2739.90719843; 

B3PW9 1/6-31 1G 

3 

-2739.93651283; 

B3PW9 1/6-31 1G 

5 

-2739.87274311; 

-2739.87274393;* 

B3PW9 1/6-31 1G** 

4 

-2740.0326636; 

B  3  PW9 1/6-31 1G 

7 

-2739.74692076; 

*  The  a  electron  from  101  orbital  is  transferee!  to  98  orbital. 


Table  5.  NMR  values  of  various  bilivcrdinc  based  quantum  computing  devices  in  various  electronic 
states  (calculated  NMR  vales  are  absolute  and  not  correlated  with  standart  TetraMethylSylane 
values). 


Method 

Bites 

Multiplic 

ity 

Atoms 

15,34 

Atoms 

16,35 

Atoms 

19,36 

Atoms 

17,37 

Atoms 

18,33 

Atoms 

20 

B3PW91/ 

6- 

311G** 

1 

1 

25.7142; 

B3PW91/ 

6- 

31 1G** 

2 

1 

25.2552; 

25.5873 

23.2951; 

B3PW91/ 

6- 

31 1G** 

3 

1 

25.9601; 

26.8833 

25.9857; 

25.6532 

23.4870; 

IJHF/6- 
31 1G** 

2 

2  (charge 
-1) 

26.8499; 

268498 

26.3916; 

B3PW91/ 

6- 

31 1G** 

4 

1 

24.7067; 

24.7482 

25.6626; 

25.5183 

25.4089; 

25.7633 

23.2992; 

B3PW91/ 
6-3 1 1 G 

6 

1 

25.6044; 

25.6044; 

25.4465; 

25.4467 

25.8578; 

25.8577 

25.5658; 

255659 

25.1355; 

251357 

26.0386; 

UHF/6- 
31 1G** 

6 

2  (charge 
-1) 

26.8151; 

26.8151 

25.8586; 

258588 

26.3127; 

263127 

25.5726; 

25.5726 

26.1080; 

261079 

26.5819; 

B3PW91/ 
6-31 1G 

6 

5 

27.7561; 

27.7566 

27.7162; 

27.7192 

28.9030; 

28.9039 

27.5793; 

27.5809 

27.4622; 

27.4644 

28.6160; 

B3PW91/ 
6-31 1G 

6 

3 

26.6736; 

26.6732 

26.2559; 

26.2558 

26.5871; 

26.5869 

26.2206; 

26.2205 

26.2687; 

26.2686 

26.7602; 

B3PW91/ 
6-31 1G 

6 

1 

26.3429; 

26.3428 

26.3639; 

26.3635 

26.5859; 

26.5864 

26.3833; 

26.3833 

26.0345; 

26.0341 

26.9946; 

B3PW91/ 
6-31 1G 

6 

5* 

27.7561; 

27.7567 

27.7161; 

27.7192 

28.9031; 

28.9038 

27.5794; 

27.5809 

27.4622; 

27.4645 

28.6160; 

B3PW91/ 

6— 

31 1G** 

6 

4** 

24.7086; 

24.0911; 

24.0908; 

24.2666; 

24.2665; 

24.1407; 

24.1408; 

24.2297; 

24.2297; 

24.5588; 

B3PW91/ 
6-31 1G 

l- . 

6 

7* 

26.0672; 

26.0680 

26.5833; 

26.5836 

26.4388; 

26.4387 

26.1883; 

26.1877 

26.3345; 

26.3348 

26.0566; 
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**  7iic  charge  is  equal  to  1  (in  calculation  was  included  Fermi  Contact  Coupling). 


Additionally  to  quantum  computing  devices  research  it  was  applied  technique  of  localized 
orbitals.  The  purpose  of  performed  work  is  obtain  place  where  unpaired  electron  is  localize. 
According  to  our  investigation  the  Cu  atom  have  not  influenced  to  shift  of  proton  NMR.  The 
chemical  shifts  of  the  investigated  molecule  depends  on  multiplicity  (see  Tables  4  and  5).  The 
increasing  multiplicity  leads  to  larger  chemical  shifts.  The  orbitals  on  which  the  unpaired  electrons 
arc  displaced  play  very  important  role.  In  the  case  of  Singlet  there  is  no  unpaired  electron  therefore 
chemical  shifts  of  the  H  atoms  are  very  similar  and  are  the  smallest  in  comparison  of  that  of  Triplet 
and  Pcnlct.  In  the  case  of  Triplet  the  N  and  C13  and  C14  atoms  possess  p,  orbitals  that  is 
perpendicular  the  plane  of  the  investigated  molecule.  On  the  orbital  the  unpaired  electron  is  located 
thus  the  chemical  shifts  of  H20  atom  that  is  closer  to  Cl  3  and  C14  is  the  laigcst  in  the  investigated 
case. 

In  the  cases  of  Pcntct  when  the  electron  from  101  orbital  is  transferred  to  98  orbital  the  above 
mentioned  p?  orbitals  arc  larger  amount  than  that  of  tiiplct.  Now  all  N  and  C8,  C9,C13,  C14,  C28, 
C27  atoms  possess  the  above  mentioned  orbital.  Thus  the  chemical  shift  of  the  H  atoms  is  the 
largest  in  comparison  of  Singlet  and  Triplet.  When  the  electron  is  not  transfer  on  from  101  to  99 
orbital,  the  electron  is  not  located  on  the  pz  orbital  of  N2  but  on  sp  type  orbital  displacement  on  the 
plane  of  the  molecule  and  directed  to  C22  atoms.  Thus,  the  chemical  shift  of  Pentets  is  coincide  nor 
the  small  changes  of  one  arc  obtained  for  HI 9  and  H36  atoms.  We  hope,  the  calculations  of  piolon 
NMR  when  it  should  be  accounted  relativistic  terms  (sec  our  developed  theory  in  the  Appendix  1) 
will  give  larger  difference. 

In  the  case  of  Septet  the  chemical  shift  of  H  is  smaller  than  that  of  Pentet  because  only  CIO, 
C27,  C29  and  C3  atoms  possess  the  pz  type  half-occupied  (one  unpaired  electron)  localize  orbital. 
Thus,  our  assuming,  that  the  chemical  shift  depends  on  displacement  of  unpaired  electron  on  C 
atoms  nearest  H  is  confirmed  from  this  localized  orbitals  technique  analysis. 

Currently  we  are  performing  programming  works  for  the  evaluation  of  relativistic  effects  in 
open  shell  species  (including  bili verdin  derivatives)  proton  NMR  calculations.  Including  correct 
interactions  of  Fermi  Contact  Coupling,  Spin— Orbit,  Spin— Spin  and  other  lelativistic  teims  (see  oui 
formulae  in  Appendix  1)  for  calculations  of  series  Cu,  Fc  or  Co  biliverdin  derivatives  and  other 
paramagnclically  'H  NMR  shifted  derivatives  should  allow  us  to  find  quantum  computer  elements 


with  broad  NMR  spectra  and  good  Quantum  Bits  resolution. 


2.3.2.  Design  of  Elements  of  Quantum  Computers  Based  on  Electronic  Paramagnetic  Resonance 


Endohcdral  fullerenc  supermolecules  ErSc2N@(CH2)(CH2)(NH)C80  possessing  from  one  to 
five  Cl  atoms  replaced  H  atoms  (sec  Figures  18-23)  were  optimized  geometries  using  ROHF 
method  and  Watanabe  basis  set.  These  six  open  shell  structures  should  possess  six  different  EPR 
values  which  should  be  possible  to  use  for  construction  of  six  QuBits  quantum  computing  elements. 


Figure  !  8.  EPR  quantum  computing  element  based  on  endohedral  lullerene  supermolecule 
ErSc2N@(CH2)(CH2)(NH)C80 


Figure  19.  EPR  quantum  computing  element  based  on  endohedral  fullerene  supermoleculc 
ErSc2N@(CHCl)(CH2)(NH)C80  possessing  one  Cl  atom  replaced  H  atom. 


Figure  20.  EPR  quantum  computing  clement  based  on  endohedral  fullerene  supermoleculc 
ErSc0N @ (CHCl)(CHCl)(NH)Cg0  possessing  two  Cl  atoms  replaced  two  H  atoms. 


Figure  21.  EPR  quantum  computing  element  based  on  endohedral  fullcrene  supermolceulc 
ErSc2N@(CHCl)(CCl2)(NH)C80  possessing  three  Cl  atoms  replaced  three  H  atoms. 


Figure  22.  EPR  quantum  computing  element  based  on  endohedral  fullcrene  supermolceulc 
ErScoN@(CHCl2)(CCl2)(NH)C80  possessing  four  Cl  atoms  replaced  four  H  atoms. 


Figure  23.  EPR  quantum  computing  element  based  on  endohcdral  fullcrcne  supcrmolcculc 
ErSc,N@(CHCU)(CCl2)(NCl)C80  possessing  five  Cl  atoms  replaced  five  H  atoms. 

Presently  we  arc  calculating  values  of  EPR  of  all  structures  presented  in  Figures  1 8-23. 

2.4.  Design  of  proton  transfer  molecular  complex  device  based  on  quantum  chemical  first 

principles  and  Ab  initio  study  of  cuprate-mediate  oxidative  P-0  coupling  of  white 

phosphorus  and  alcohols 

Several  years  ago  Ya.  Dorfman  et  al  [11,  12]  created  new  technology  based  on  selective  and 
fast  reaction  in  which  white  phosphorus  react  directly  with  organic  substrates.  In  the  above- 
mentioned  technology  the  catalytic  activation  and  functionalization  of  white  phosphorus  in 
alcohol/arcnc  solutions  of  Cu(Il)  and  Fc(III)  salts  under  mild  reaction  conditions  and  dioxygen  as  a 
cheap  and  efficient  oxidant  arc  used.  This  reaction  mechanism  should  be  useful  for  design  of  the 
proton  transfer  molecular  devices.  Therefore,  a  variety  of  Cu  complexes,  that  are  thought  to  be 
relevant  models  for  representing  proton  transfer  complexes  have  been  investigated  by  first 
principles  and  ab  initio  quantum  chemical  methods.  We  assumed  that  in  our  investigated  complexes 
the  alcohol  molecule  transfer  proton  what  is  accepted  by  one  ligand  of  the  coordinated  complexes. 
The  phosphorus  molecule  could  play  role  of  acceptor  of  deprotonated  alcohol.  Thus,  if  the  pioton 
take  place  on  alcohol,  one  could  have  logical  0  (or  switch  out),  while  the  remove  of  proton  from 
alcohol  indicate  logical  1  (or  switch  on). 

Firstly  the  coordinated  complexes:  1)[CuC12NH3CH3OH  n2-P4]°,  2)|CuC1CH3OH(NH3)3  n2- 
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P4]+l  and  3)  [CuCH3OH  (NH3)3  n2-P4]+2  that  are  relevant  in  the  Ya.  Dorfman  et  al.  created 
technology  were  formed  and  investigated.  The  alcohol  is  involved  in  the  coordinated  complexes 
due  to  assumption  of  intramolecular  proton  transfer  i.e.  Only  the  alcohol  molecule  belonging 
coordinated  sphere  could  lost  proton.  The  cuprate-ligand  bond  lengths  have  been  considered 
approximately  equal  to  the  sum  of  the  covalent  radii  [13].  The  initial  starting  geometries  of  the  1-3 
complexes  were  adjusted  applied  Density  Functional  Theory  Method  B3PW91  [6]  in  the  6—31 1G" 
basis  set  with  full  Bcrny  geometry  optimization  until  the  stationary  points  on  the  potential  surfaces 
were  found.  Then  the  Hessian  of  the  above  complexes  is  investigated  to  be  sure  that  the  real 
stationary  point  is  obtained.  The  proton  transfer  devices  (switches)  is  designed  basing  on  the  results 
of  investigation  the  overlap  populations  (OP),  total  atomic  charge  and  localize  orbital  investigation. 
A  positive  value  of  OP  points  to  the  presence  of  a  bonding  between  atoms,  while  a  negative  OP 
value  indicates  antiboding  that  are  confirmed  by  analysis  of  bond  orbital.  The  bond  orbital  analysis 
allows  us  to  obtain  atoms  that  arc  capable  to  form  new  chemical  bonds. 

As  it  was  mentioned  above  the  goal  of  our  study  is  exhibit  that  the  proton  could  be  tiansfencd 
at  the  medium  of  transition  metal  complexes  that  arc  efficient  oxidant.  Firstly  the  Berny  geometry 
optimization  procedure  applying  B3PW91/6-311G**  method  of  the  1,  2  and  3  complexes  was 
performed.  According  to  the  investigation  the  Cu  metal  with  the  ligand  forming  square-planar 
complexes  (Figure  24).  The  angles  and  distances  between  the  atoms  that  form  the  square-planar 
complexcs  arc  presented  in  Table  6. 


Figure  24.  The  geometry  of  the  1  (first)  complex. 


The  distances  between  the  phosphor  atoms  are  larger  than  the  sum  of  covalent  tadii.  On  the 
other  hand,  the  OP  values  between  the  above  atoms  is  approximately  equal  to  zero  that  indicate 

absence  of  any  bonding.  It  implies,  that  the  P4  molecule  does  not  belong  the  coordination 
complexes  of  Cu  metal.  Thus  assumption  of  the  intramolecular  proton  transfer  is  not  conventional. 
In  keeping  this  results,  we  must  have  in  mind  possibility  of  the  intermolecular  proton  transfer. 

Table  6a.  The  distances  (A)  between  atoms  in  the  square-planar  complexes. 


Atoms/Distances 

1  complex 

2  complex 

3  complex 

Cu-Cl  (N) 

2.22 

2.2 

2.01 

Cu-0 

2.09 

2.03 

2.05 

Cu-Cl  (N) 

2.22 

2.03 

2.01 

|  Cu-N 

2.01 

2.02 

2.01 

Table  6b.  The  angle  (degrees)  between  bonds  of  the  coordinated  complexes.  The  numbcis  ol 
atoms  is  presented  in  Figures  24-26. 


Bonds  of  complex 

1 

Angles 

Bonds  of 

complex  2 

Angles 

Bonds  of  complex 

3 

Angles 

07-Cul-N14 

88.09 

06-Cul-C113 

85.07 

06-Cul-N12 

87.5 

07-Cul-C113 

82.76 

06-Cul-C112 

90.63 

O6-Cul-N20 

87.67 

Cl  1 6— Cu  1  —Cl  1 3 

101.74 

C112~Cul-N17 

88.2 

N16-Cul-N12 

92.09 

Cl  1 6-Cul-Cl  14 

88.09 

N13-Cul-N17 

96.1 

N16-Cul-N20 

92.46 
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Figure  26.  Geometry  of  the  complex  3. 


Figure  27.  Proposed  reaction  mechanism. 


Theoretical  investigations  proved  that  some  P-P  bonds  are  very  weak  or  destructed  in 
comparison  to  that  of  isolated  tetraphosphorus  molecule  (see  Table  7).  The  investigation  ol  localize 
orbital  of  the  above  atoms  indicate  presence  of  sp  type  orbitals. 


Table  7.  The  overlap  population  values  between  P  atoms  in  the  isolated  tetraphosphorus 
molecule  and  one  in  this  molecule  near  coordinated  complexes  Cu. 


Bonds 

P4 

Complex  1 

Complex  2 

Complex  3 

P2-P3 

0.13 

0.03 

0.08 

-0.14 

P2-P4 

0.13 

0.09 

-0.01 

-0.15 

P2-P5 

0.13 

0.06 

0.01 

-0.14 

P3-P4 

0.13 

0.16 

0.17 

0.26 

P3-P5 

0.13 

0.15 

0.18 

0.26 

P4-P5 

0.13 

0.16 

0.18 

0.26 

It  implies,  that  the  single  tetrahedral  phosphor  molecule  is  decomposed  to  P  atom  and  P3  ion 
at  the  medium  of  this  coordinated  complexes  could  form  new  chemical  bond  and  it  is  very 

important  for  successful  nucleophilic  attack  of  the  P4  molecule  by  the  alkoxide  ion.  Moreover,  the 

total  atomic  charge  on  some  atoms  in  P3  is  positive  that  is  a  favorable  condition  to  account  for  the 
P-0  coupling  also  (Table  8). 

Thus  the  destructed  P4  molecule  behaves  as  a  strong  acceptor  of  negative  charge  and,  at  least 


in  part,  this  feature  explains  how  one  of  the  phosphor  atom  efficiently  undergoes  the  attack  of  the 
methoxide  ion. 


Table  8.  The  charge  on  the  atoms  of  the  isolated  tetraphosphorus  molecule  and  one  on  the 
molecule  near  coordinated  complexes. 


Atoms 

Complex  1 

Complex  2 

Complex  3 

P2 

-0.13 

-0.16 

-0.43 

P3 

0,1 

0.13 

0.13 

P4 

-0.03 

-0.01 

0.14 

P5 

-0.03 

0.01 

0.12 

j  P, ions 

0.06 

0.21 

0.39 

The  Cl  atoms  and  NH3  group  are  strongly  coordinated  to  Cu  in  all  described  complexes  (OP 
value  ranging  from  0.1 1  to  0.29).  A  weak  intermolecular  interaction  is  established  between  Cl  and 
the  proton  HI 2  or  HI  1  of  alcohol  in  the  1  and  2  complexes  respectively.  The  OP  value  is  equal  to 
0.08  and  0.02.  The  locali/.e  orbital  analysis  indicate  that  in  the  1  complexes  Cu  and  ligands  made 
coordinational  bond.  The  Cl  13  atoms  possess  the  partly  filled  (one  electron)  orbital  directed  to 
PI  12  atom.  The  other  Cl  atom  of  the  1  complexes  have  the  above  mentioned  orbital  also  but  it  is 
directed  out  of  coordinated  sphere  of  the  transition  metal.  A  non-innocent  role  can  therefore  be 
assigned  to  the  halide  ligand  which  behaves  as  weak  bases  capable  of  accepting  protons  from  the 
alcohol.  Such  hydrogen-bonding  interaction  between  the  Cu-coordinated  chloride  and  the  proton 
represents  a  favorable  condition  to  account  for  the  deprotonation  of  the  alcohol  ligand  and  the 
transfer  to  alkoxidc  ion  to  the  proximal  phosphorus  atom.  The  ultimate  result  ol  this  likely 
concerted  process  is  the  formation  of  P-0  bond  that  is  not  realized  in  the  3  complex  due  to  absence 
halide  and  as  consequence  docs  not  realize  the  alcohol  deprotanation.  It  implies,  that  halide  ligand 
do  not  directly  participated  in  the  oxidative  P-0  coupling  but  they  are  active  components  of  the 
reaction  providing  nucleophilic  assistance  to  the  alcohol  deprotonation  step  as  it  was  obtained  in 
i  141. 

Notwithstanding  favorable  requisites  to  promote  the  deprotonation  step  of  the  coordinated 
alcohol,  the  overall  simulation  on  the  Cu  compounds  1  and  2  point  to  a  stronger  coordination 
between  the  Cu  and  CH3  anion.  The  OP  values  of  Cu-0  in  complex  1  is  equal  to  0.08  and  it 
becomes  large  for  complex  2  (0.12).  The  localize  orbital  analyses  indicate  strong  coordination  also. 
These  strong  coordination  bonds  make  unlikely  the  occurrence  of  an  intramolecular  nucleophilic 
attack  of  the  alkoxidc  group.  In  keeping  with  this  results,  it  is  possible  to  foresee  that  pioton 
transfer  could  not  realized  in  our  designed  2  complexes  while  it  could  be  present  in  the  1  complexes 
as  intermolecular  due.  to  presence  of  two  halide  capabling  deprotonate  two  alcohols. 


As  indicated  our  latest  results  of  geometry  optimization  and  Hessian  investigations,  the 
[CuCl  NH’,CH;,0H1()  square-planar  complex  could  be  forming  or  existing  without  presence  ol 
tetraphosphorus.  The  obtained  geometry  coincides  with  that  of  the  1  complexes.  The  results 

confirm  our  conclusion  that  P4  do  not  belong  for  coordinated  sphere  of  Cu.  In  keeping  with  the 
above  results,  we  proposed  that  in  the  copper-mediate  the  tetraphosphorus  molecule  is  destructed 
and  two  halide  belonging  the  coordinated  sphere  of  the  transition  metal  atom  could  be  able  to 
dcprolonalcd  two  alcohol  molecule  one  of  them  is  not  belonged  the  coordinated  sphcic.  This 
assuming  is  confirmed  by  bond  orbital  analysis.  According  to  our  investigations  the  Cl  atom,  that 
could  capable  to  deprotonate  the  non-coordinated  alcohol  molecule,  possess  four  orbitals.  One  of 
them  indicate  non  covalent  bond  between  the  Cl  atom  and  Cu  (therefore  OP  value  between  these 
atoms  is  small  positive),  other  is  composed  by  the  s-type  orbital  (78%)  and  p-type  orbital  (22%), 
and.  lastly,  two  p-type  orbitals  placement  in  the  x-y  and  z-y  planes  and  directed  from  coordinated 
sphere  of  Cu.  The  last  mentioned  orbitals  could  form  new  chemical  bonds.  Thus,  this  non- 
coordinated  anion  make  the  occurrence  of  an  intermolecular  nucleophilic  attack  of  the  alkoxide 
group  in  the  ease  of  complex  1  and  support  the  experimental  results  suggesting  that  the 
phosphorylation  reaction  proceeds  only  in  the  presence  of  Cu(II)  to  give  alkylphosphates. 

Finally,  all  our  observations  indicate  that  the  1  complex  is  nicely  tailored  for  modeling  the 
proton  transfer  devices.  Summarizing  all  results  we  may  proposed  mechanism  of  the  above  reaction 
(Figure  27):  in  the  Cu  medium  this  metal  square-plane  coordinated  complexes  is  formed  and 
destructed  the  tetraphosphorus  molecule;  halide  ligand  do  not  directly  participate  in  the  oxidative 
P-0  coupling  but  they  arc  active  components  of  the  reaction  providing  nucleophilic  assistance  to 
the  alcohol  (that  docs  not  depend  the  coordinated  sphere  of  Cu)  deprotonation  stcp;thc 
intermolecular  nucleophilic  attack  of  alkoxide  group  is  piesent. 

In  keeping  these  results  of  investigation  the  complexes  4)  [CuCMNHj^l0  and  5) 
[CuCF(NH3)2l0  were  designed  as  more  favorable  for  proton  transfer  devices.  In  the  5  complex  one 
Cl  atom  possesses  p  type  orbital  that  capable  to  form  new  chemical  bond  i.  e.  could  accept  proton 
from  alcohol.  But  in  the  case  the  acceptor  of  deprotonated  alcohol  molecule  absent,  thus  switching 
devices  could  be  obtained  as  increasing  concentration  of  CH30-  or  polarity  of  media. 

CONCLUSIONS: 

1.  OR  logically  controlled  molecule  machine  designed  from  benzene-N=N-benzene-N02  with 
electron  acceptor  fragment  -(PI1-NO2)  (moving  part  after  excitation)  and  two  photoelectron  donor 
pans:  dithicno|3,2-b:2A3’-d]lhiophene  and  ferocene  molecules  show  some  not  intensive 
absorption  in  infrared  region.  Applied  DFT-TD  method  and  our  modified  visualization  piogiam 
showed  small  charge  transfer  in  first  and  second  (IR)  excited  states. 


2.  OR  logically  controlled  molecule  machine  designed  from  thiobenzene-N=N-benzene-N02 
with  electron  acceptor  fragment  -(Ph-N02)  (moving  part  after  excitation)  and  two  photoelectron 
donor  parts:  dilhicno|3,2-b:2’,3’-dlthiophene  and  ferocene  molecules  was  optimized  by  HF 
method  and  shows  characteristic  (similar  to  tetrahedral)  valence  angles  near  the  S  atom  inserted  in 
the  six  membered  ring. 

4.  It  was  designed  three  variable  logically  controlled  molecule  machine  for  the  photoinduced 
moving  in  3Dimcnsional  zeolite  surfaces  and  construction  of  3D  addressing  memory. 

5.  Design  and  optimization  of  geometry  by  HF/6-31G  was  performed  lor  two  different  selective 

two  variable  molecular  logic  machines. 

6.  It  was  designed  two  different  binary  two  and  three  variable  molecular  logic  machines  moving 
in  accordance  to  relaxation  energies  due  to  various  photo  excitations. 

7.  We  performed  design  and  calculations  of  three  variable  molecular  logical  devices  based  on 
organic  electron  donor:  dithicno|3,2-b:2\3’-d]thiophene,  TTF  and  ferocene  (CioHioFe)  and 
electron  acceptor  molecules:  1 ,3-bis(dicyanomcthylidene)indanc  and  endohedral  fullerene 
ErSc2N@C80  substituted  derivative  ErSc2N@(CH2CH2NH)C80,  electron  donor-bridge-electron 
acceptor  dyads  and  triads  including  electron  donor  and  acceptor  molecules  joined  with  -CH=CH- 
bridec.  Design  of  new  series  molecular  implementations  of  two  variable  logic  functions:  AND 
(NAND).  OR  (NOR)  is  based  on  quantum  chemical  ab  initio  HF/6-311G  geometry  optimization 
procedure  which  shows  that  our  newly  designed  logical  gates  based  on  ErSc2N@C80  substituted 
derivative  ErSc2N@(CH2CH2NH)C8()  are  more  stable  rather  than  based  on  ErSc2N@C80 
substituted  derivative  ErSc2N@(CH2)C80- 

8.  It  was  performed  design  and  calculation  of  biliverdine  derivatives  which  possess  from  one  to 
six  Quantum  Bits  generating  by  proton  NMR  for  closed  and  open  shell  species  (including  Fcimi 
Contact  interaction  but  not  including  other  relativistic  effects).  Theory  of  our  obtained  diffcient 
relativistic  terms  for  open  shell  paramagnctically  shifted  proton  NMR  calculations  and  formulae 
prepared  for  programming  and  implementation  to  the  DALTON  package  arc  given  in  Appendix  1 
of  this  Report  04. 

9.  Set  consisting  from  six  various  ErSc2N@(CH2CH2NH)C80  derivatives  where  H  is  replaced 
by  Cl:  ErSc2N@(CH2CH2NCl)C80,  ErSc2N@(CH2CClHNH)C80,  ErSc2N@(CClHCH2NH)C80, 
ErSc2N@(CH2CCl2NH)C80  etc.  were  designed  and  optimized  by  FIFAVatanabe  in  order  to  get 

quantum  computing  elements  possessing  six  QuBits  based  on  six  different  EPR  signals. 

10.  It  was  designed  molecular  switch  based  on  proton  transfer  in  molecular  complex  device 
composed  from  cuprate-mediate  oxidative  P-0  coupling  of  white  phosphorus  and  alcohols. 
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APPENDIX  1 

Theoretically  Obtained  Relativistic  Terms  for  Proton  NMR  Calculations  of  Open  Shell 
Paramagnetic  Species  rewritten  in  more  clear  and  exact  form  and  prepared  for  the  programming 
works  for  the  implementation  in  DALTON  package: 


n 


I.  Electron  Zeeman  operator. 


«* = £s(.>B-£feA + XA  +5.A  )=  I>,A + &A +X>A  =  SA + s^, + s,s. 


II.  Fermi-Contact  operator. 


H iC  =  £S(0  ■  /S (r,  -  r„)  =  (?7,  +  S,I,  +  S7, £6(1;  -  r„)  =  ($7,  +  Sfl,  +  57, 


III.  Response  function  «FC;SZ». 
I  will  write  down  matrix  elements: 
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I  used  selection  rules  from  ncwelemcts.doc  for  matrix  elements  to  remove  equal  to  zero 
elements.  In  same  way  I  will  get: 


((|/ir,,jn}(4«,4d  +s,/,  +s:/i  H‘-S*-Sv +5cal()  ={(K/a4"K"1s^I(I +((^z,5ri4"K'1KB'l4 

={45zfivJ»X»MK} + 


Now  1  will  define  response  function  «FC;SZ»  energy  terms  ESz-fc: 


where  clement  Eab-  correspond  to  energy  term  with  a  and  b  cartesian  components  of  B 
and  I  respectively.  In  same  way  I  will  get  tern  corresponding  to  Efc-sz: 
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Now  I  will  derive  nuclear  shield  tensor  components  for  «SZ;FC»  response  function: 
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IV.  Spin-Dipolar  operator. 
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V.  Response  function  «SD;SZ». 
I  will  write  down  matrix  elements: 
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In  same  way  I  will  gel: 
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Now  1  will  define  response  function  «SD1;SZ»  energy  terms  Esz-sdo 
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In  same  way  I  will  gel  tem  corresponding  to  ESm-sz: 
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For  Hsd2  I  will  obtain  results  in  same  way  as  for  Hrc: 
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Now  I  will  derive  nuclear  shield  tensor  components  for  «SZ;SD2»  response  function: 
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